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Field theory calculations predict multiplicative logarithmic corrections to correlation functions from 
marginally irrelevant operators. However, for the numerically most suitable model - the spin-1/2 
chain - these corrections have been controversial. In this paper, the spin-spin correlation function 
of the antiferromagnetic spin-1/2 chain is calculated numerically in the presence of a next nearest 
neighbor coupling J2 for chains of up to 32 sites. By varying the coupling strength J2 we can 
control the effect of the marginal operator, and our results unambiguously confirm the field theory 
predictions. The critical value at which the marginal operator vanishes has been determined to be 
at J 2 /J = 0.241167 ± 0.000005. 
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The spin-1/2 chain has attracted much attention in 
theoretical physics since the early days of quantum me- 
chanics as a simple model for many-body effects. How- 
ever, correlation functions in the antiferromagnetic chain 
could only be predicted with the help of field theory 
calculations [g], and it was only a few years ago when 
a prediction of a multiplicative logarithmic correction 
to the asymptotic spin-spin correlation functions was 
made coming from a marginally irrelevant operator in 
the field theory description ||]4| . At first sight this result 
seems surprising, since it indicates that the effect of a 
marginally irrelevant operator effectively increases in the 
long distance limit, while on the other hand its effect on 
other quantities (energy spectrum [gj, susceptibility j^]) 
vanishes as length L — > 00 and temperature T — > 0. 

However, the underlying calculation which predicted 
the multiplicative correction from a marginal operator is 
quite general || and of great importance in other mod- 
els as well (e.g. the 4 -model). Since the spin-1/2 chain 
is accessible by a number of numerical methods it rep- 
resents an ideal testing model for this important result. 
Therefore, considerable effort has been made to verify 
the logarithmic corrections (7j|||^, but often different or 
no multiplicative corrections were found. 

In the present paper we choose to numerically diago- 
nalize the model Hamiltonian of up to 32 sites exactly. 
While this method is limited by relatively small system 
sizes, it allows us to introduce a next nearest neighbor 
coupling J 2 quite easily and, moreover, the results are 
very accurate. It is known that by adjusting J2 we can 
change the bare coupling strength of the marginal opera- 
tor and even set it to zero without destroying the validity 
of the field theory continuum limit [jlO| . Since we are able 
to probe the system with a reduced marginal coupling, 
the range of validity for the perturbative field theory cal- 
culations is increased exponentially. Moreover, we will 



use a novel approach to analyze the data which does not 
depend on any extrapolation to an infinite system size 
for the field theory calculations to be valid, so that 32 
sites turn out to be very sufficient to confirm the effect 
unambiguously. 

Our model Hamiltonian is the antiferromagnetic spin- 
1/2 chain (J > 0) with a next nearest neighbor coupling 
■h 



h — j s x ■ s x+ i + J2 S x ■ S 



x+2- 



(1) 



In the long- wavelength, low-energy limit this model is 
well understood in terms of the Wess-Zumino-Witten 
(WZW) non-linear a-model with a topological coupling 
constant k= 1 pL In the continuum limit, the spin-1/2 
operator at each lattice site x can then be expressed 
in terms of a WZW SU(2) matrix and (related) chiral 
SU(2) currents 3l,r 



const. i(—l) x tr[gcr]. 



(2) 



Strictly speaking this theory is only valid up to cor- 
recting higher order operators in the Hamiltonian which 
we have neglected so far. In particular, we want to con- 
sider the effect of the leading irrelevant marginal operator 
which can be written in terms of the current operators 
with some coupling constant A 



SH 



-2-kXJl ■ Jr. 



(3) 



For periodic boundary condition the next "leading" irrel- 
evant operator in the Hamiltonian has scaling dimension 
4, which can be neglected even for modest lengths of or- 
der L > 10 and temperatures of order T < 0.1J JT^ |. 

The effect of the marginal operator on the uniform (i.e. 
current-current) part of the spin-spin correlation function 
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is well understood in terms of a simple additive logarith- 
mic correction. This leads to interesting predictions for 
the asymptotic behavior of the susceptibility which 
were well confirmed by Bethe ansatz and also experimen- 
tal results The situation is quite different, however, 
for the alternating part of the spin-spin correlation func- 
tion G(r) which is given in terms of the WZW field g 



G(r) oc< tr[ ff er](0) • tr[ga](r) > 



(4) 



For the primary field g we have to take into account the 
effect of the marginal operator on the anomalous scal- 
ing dimension 7(A) of tr[ger] as well (while the scaling 
dimension of the currents is fixed) . The correlation func- 
tion then also depends on the marginal coupling constant 
A, and G(r, A) obeys the renormalization equation 



-<L +/3 (A)£ + 27(A) 



G(r, A) = 



(5) 



where /3(A) is the beta function of the marginal operator 
STL. This equation (^]) is then solved by introducing a 
scale dependent coupling constant A(r) which obeys 
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/3(A) 



(6) 



and writing (up to a prefactor which may contain additive 
corrections in A) 



G(r) oc exp 



exp 



In T 

2 / j(x)d(lnx) 

/ In n 

"AM 7(A) 



An 



/3(A) 



dX 



(7) 



where r% is an ultraviolet cut-off of the order of the lat- 
tice spacing and Ao is the unrenormalized (bare) coupling 
strength at that energy-scale. The beta-function and the 
anomalous dimension are known as perturbative expan- 
sions in A IH1 



/3(A) = -A 2 + 0(A 3 ) 

iW = \-\ + o{x 2 ) 



(8) 



Integrating equation (|6j), we obtain an expression for A 
up to 0[ln(lnr)/(lnr) 



A(r) 



An 



l + A ln(r/n) ln(r/r )' 



(9) 



and together with equation (]?]) we can determine the cor- 
relation function up to higher orders in A 

G(r) oc (A/A ° rl/2 [1 + 0(A)] « ^ A ° lnr/r ° (10) 
r r 

From equation (^) we see that the constant rn is simply 
given by 



rn = nexp(— 1/Ao), 



(11) 



In the asymptotic limit we find G oc Vhi r/r as al- 
ready predicted in reference || unless Ao is infinitesimally 
small. Therefore, very close to the critical point the pre- 
dicted asymptotic limit will never be reached although 
the perturbative expansions in equation (|^) are very ac- 
curate. Since rn is a non-universal constant, it is essential 
to keep the full expression in equation ( |To| ) to test that 
result (as opposed to just assuming the asymptotic limit). 

It is important to note that equation (^J) is only valid 
if r is the only length scale in the system, but all previ- 
ous studies assumed that this would require the infinite 
length limit L — ► 00 which is numerically very difficult to 
achieve. However, it is perfectly feasible to study the sys- 
tem at finite, but varying lengths L and keeping the ratio 
R = r/L fixed. Since r = RL is no longer independent 
of L, there is only one variable parameter in the problem 
and therefore expression © holds for G(r = RL) (for 
each R separately, i.e. the overall proportionality con- 
stant will depend on the choice of R). To test the field 
theory predictions it is therefore much easier and more 
accurate to select the fixed ratio to be R = 1/2, and in 
what follows we will hence set 



L/2. 



(12) 



It is then straightforward to determine the correlation 
function G(r — L/2) numerically for different lengths 
L < 32 and coupling constants J2 (as opposed to se- 
lecting R = and having to choose some controversial 
extrapolation to the L — > 00 limit as in all previous 
studies). 

The correlation functions of the ground state were 
found by exact diagonalization using the modified Lanc- 
zos method on periodic spin chains of up to 32 sites (using 
ca. 4.7 million basis states). Table | shows the results for 
the spin-spin correlation functions < S z (0)S z (r) > for 
different distances r — L/2 and for various values of the 
next nearest neighbor coupling J2. 

By adjusting J 2 we can change the bare coupling An 
and even set it to zero at some critical value J" 1 ' [^0| . 
To lowest order, Ao increases linearly with the differ- 
ence from the critical point A J2 = J™' — J2 ■ The critical 
value J™* was determined to be 



J. 



crit 



0.241167 ±0.000005 J, 



(13) 



which we believe to be the most accurate estimate to 
date. This result was obtained by examining the energy 
difference A between the first two excited states of total 
spin S = and S = 1 in the periodic chain spectrum 
as a function of length L and J2. Those two states are 
nearly degenerate and their energy difference is only due 
to the marginal operator and higher order terms in the 
field theory Hamiltonian. 
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r = L/2 


,J 2 = -0.25J 


J 2 = 


J 2 = 0.1 J 


J 2 = 0.2J 


J 2 = J 2 crit 


4 


0.0572167 


0.0497077 


0.0459760 


0.0417541 


0.0398672 


5 


-0.0545609 


-0.0469513 


-0.0433261 


-0.0390852 


-0.0370562 


6 


0.0423313 


0.0356626 


0.0325190 


0.0288716 


0.0271334 


7 


-0.0403820 


-0.0340938 


-0.0311410 


-0.0276488 


-0.0259367 


8 


0.0336161 


0.0279325 


0.0252780 


0.0221537 


0.0206302 


9 


-0.0322537 


-0.0269130 


-0.0244137 


-0.0214326 


-0.0199523 


10 


0.0279427 


0.0230233 


0.0207281 


0.0179995 


0.0166504 


11 


-0.0269492 


-0.0223039 


-0.0201300 


-0.0175197 


-0.0162116 


12 


0.0239555 


0.0196206 


0.0175960 


0.0151709 


0.0139598 


13 


-0.0232004 


-0.0190837 


-0.0171551 


-0.0148271 


-0.0136519 


14 


0.0209966 


0.0171185 


0.0153044 


0.0131186 


0.0120183 


15 


-0.0204030 


-0.0167013 


-0.0149649 


-0.0128595 


-0.0117902 


16 


0.0187105 


0.0151987 


0.0135530 


0.0115605 


0.0105511 


TABLE I. The spin-spin correlation function 


< S z (0)S z (r) > for 


various values of r 


= L/2 and J 2 . 


In some cases the last 



digit is uncertain. 



Wc find the critical point by determining the value of 
Ji at which the energy difference A is exactly propor- 
tional to 1/L 3 , because at that point the correction A 
comes only from the higher order terms and therefore the 
marginal operator is absent [ jl2| . This approach proves 
to be more accurate than the conventional method of ex- 
trapolating the value of J2(A = 0) as L — > oo (which of 
course gives the same value, although with less accuracy 
@). ' 

As we approach J| rit the value of the effective length 
scale ro becomes exponentially small ro = r\ exp(— 1/Ao), 
and the expansions in l/ln(r/ro) in equations (^|) and 
( |io| ) give very accurate results even for moderate lengths. 
However, we require the marginal operator to remain 
the leading correction (compared to the 1 / L 2 corrections 
from higher order operators), so it is not always useful to 
examine the system arbitrarily close to J| rlt . 
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a J 2 =-0.25J 
i J,=0 
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FIG. 1. The multiplicative correction rG(r) as a function 
of r for various values of J 2 according to the numerical data 
in table I. 
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FIG. 2. The fitting curves to rG(r) according to equation 
([uj) with the parameters given in table [o]. 

To extract the multiplicative correction we perform a 
cubic spline interpolation of r < S z (0)S z (r) > from the 
data in table | for even and odd r separately. The differ- 
ence between the two curves then gives the multiplicative 
correction rG(r) to the alternating part of the correlation 
function as shown in figure (Q) for different values of J2. 

Figure (||) shows the excellent asymptotic fit to the 
predicted form of the multiplicative corrections 



rG(r) = \J a\ \nr/r 



(14) 



with the fitting parameters given in table |H| where a 
is a normalization constant which must be independent 
of Ao. In all cases the data approaches the asymp- 
totic curve very quickly. The small deviations are pos- 
itive for small values of AJ2 = J™' — ^2 and neg- 
ative for larger AJ2, which indicates a partial can- 
cellation of higher order terms with varying relative 
magnitude as Aq is changed (this explains why the 
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first order calculation already gives such good results). 



J2/J 


In r 


a\ 


-0.25 


-1.4976 


0.021601 


-0.1 


-2.1683 


0.014864 





-2.9501 


0.010696 


0.05 


-3.6111 


0.008626 


0.1 


-4.6883 


0.006546 


0.15 


-6.9080 


0.004383 


0.2 


-14.393 


0.002080 



TABLE II. The fitting parameters for according to equa- 
tion ([li]) for selected values of J2. 



For the special J 2 = case we find general agreement 
with reference |t]] where a similar fit was made, however 
with much larger deviations, and - more importantly - 
their data does not approach the predicted form asymp- 
totically as it must. (In that reference an extrapolation 
scheme with additional adjustable parameters was used 
as well as a less accurate numerical method.) 








FIG. 3. The logarithm lnro vs. 1/aAo from table |l| shows 
an excellent fit to lnro = lnri — 1/Ao with n = 0.85 and 
a = 0.0296, confirming equation jl^). 

As the critical point is approached, we require simul- 
taneously ro — > and Ao — > in such a way that the 
correlation function G(r) remains finite. In particular 
we know from equation ( pT| ) that r and Ao are related 
(to lowest order the correlation function G(r) is indepen- 
dent of the size of the marginal coupling at the ultraviolet 
cutoff ri) i.e. for all J2 < J| rlt we must have: 

A lnn/ro = l (15) 

As can be seen in figure (|3|) this relation holds very ac- 
curately for the parameters in table [0] for all values of J 2 
with the constant in equation ( fil| ) given by a = 0.0296. 
Equation fll5|) also fixes the ultraviolet cutoff n = 0.85 



at which the bare coupling constant is defined. Because 
a and r\ are fixed, Ao and ro are related by equation JlE|), 
and there remains really only one free parameter for each 
of the curves in figure (|^) , which makes the quality of the 
fits even more impressive. 



0.8 




Moreover, analyzing the fitting parameters in table O, 
we can determine the functional behavior of the bare cou- 
pling Ao in form of a series expansion in the microscopic 
parameter A J 2 — J™* — J 2 

A = C1AJ2 + c 2 Ajf + c 3 Ajf , (16) 

As shown in figure this relation ( fL6| ) seems to 

hold even for relatively large values of AJ with c\ ~ 
1.723/J, c 2 w -1.35/J 2 , and c 3 « 1.76/J 3 . Therefore, 
once the constants a, n, and Ci are known, the shape 
and overall magnitude of G(r) can be predicted accu- 
rately for a large range of values J2 and r > 5, which is 
a clear success of the field theory calculation. 

In conclusion, we have unambiguously shown that the 
field theory calculations for multiplicative corrections to 
the correlation functions in the spin-1/2 chain are valid. 
The constants a, r\, and Ci were determined, and we 
were able to improve the estimate of J" 1 *- The basic 
result carries over to other models with marginal opera- 
tors (1-D Hubbard model, other half- integer spin chains, 
4 -models) where we also expect that the correspond- 
ing field theory calculation gives accurate predictions to 
multiplicative corrections to correlation functions. 
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